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ABSTRACT:  This  work  investigates  the  fundamental  nature  of  sulfur  mustard  surface  •  Cl  •  C  •  H 

adsorption  by  characterizing  interfacial  hydrogen  bonding  and  other  intermolecular 
forces  for  the  surrogate  molecule  (simulant)  2-chloroethyl  ethyl  sulfide  (2-CEES). 

Adsorption  at  the  surface  of  amorphous  silica  is  the  focus  of  this  work  because  of  silica’s 
low  chemical  reactivity,  well-known  properties,  and  abundance  in  the  environment.  2- 
CEES  has  two  polar  functional  groups,  the  chloro  and  thioether  moieties,  available  to 
accept  hydrogen  bonds  from  free  surface  silanol  groups.  Diethyl  sulfide  and 
chlorobutane  are  also  investigated  to  independently  assess  the  role  of  the  chloro  and 
thioester  functionalities  in  the  overall  adsorption  mechanism  and  to  explore  the 
interplay  between  the  charge  transfer  and  electrostatic  contributions  to  total  hydrogen- 
bond  strength.  Our  approach  utilizes  infrared  spectroscopy  to  study  specific  surface- 
molecule  interactions  and  temperature-programmed  desorption  to  measure  the 
activation  energy  for  desorption  of  adsorbed  molecules.  Our  results  indicate  that  2-CEES  adsorbs  to  silica  by  hydrogen 
bonding  through  either  the  chloro  or  thioether  moieties  but  is  unable  to  form  a  more  stable  configuration  in  which  both  polar 
groups  interact  simultaneously  with  adjacent  silanol  groups.  The  activation  energy  for  desorption  of  2-CEES  is  nearly  43  kj/mol, 
driven  by  both  strong  hydrogen  bonding  and  other  non-bonding  interactions.  A  systematic  study  of  chloroalkanes  reveals  that 
each  methylene  group  contributes  approximately  5—8  kj/mol  to  the  overall  desorption  energy. 


1.  INTRODUCTION 

Sulfur  mustard  [bis(2-chloroethyl)  sulfide,  HD]  is  an  extremely 
toxic  compound  that  saw  widespread  use  in  the  First  World 
War  as  a  vesicant  chemical  warfare  agent  (CWA).1  Although 
the  Chemical  Weapons  Convention  (CWC)  in  1997  mandated 
the  destruction  of  mustard  stockpiles  worldwide,  HD  can  be 
readily  synthesized  by  terrorists  or  militaries  who  choose  to 
ignore  the  CWC  policies.2  Furthermore,  HD  is  known  to  be 
persistent  in  the  environment  and  can  remain  at  sufficiently 
high  concentrations  to  |>ose  serious  health  risks  several  years 
after  initial  deployment/  Motivated  by  this  threat,  researchers 
have  directed  significant  effort  toward  the  design  of  effective 
protection  strategies  including  sorbent  materials,  decontamina¬ 
tion  procedures,  and  sensors.  These  technologies  often  rely  on 
molecule— surface  interactions  or  reactions  that  are  tailored  to 
achieve  a  particular  outcome.  However,  surprisingly  little  is 
known  about  the  fundamental  nature  of  the  surface  chemistry 
of  sulfur  mustard. 

As  HD,  a  moderate  vapor  pressure  liquid  at  room 
temperature,  interacts  with  the  surface  of  a  material,  it  may 
accommodate  to  the  surface  prior  to  desorption,  diffuse  into 
the  bulk  of  a  porous  material,  or  undergo  a  chemical 
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transformation.  The  pathways  that  govern  a  particular  HD— 
surface  interaction  are  governed  in  large  part  by  the  initial 
binding  energy  of  the  molecule,  which  controls  the  molecular 
residence  time  on  the  surface.  For  molecules  like  HD  that  have 
multiple  functional  groups  with  available  lone  pairs  of  electrons, 
one  of  the  most  likely  initial  binding  pathways  is  through  the 
formation  of  hydrogen  bonds  with  surface  hydroxyl  groups.  For 
example,  interfacial  hydrogen  bonding  has  been  shown  to  be 
the  critical  first  step  in  the  decomposition  of  2-chloroethyl  ethyl 
sulfide  (2-CEES  or  “half-mustard”)  on  the  surface  of  a 
hydroxylated  Si02— Ti02  composite  nanomaterial.4’5  Further¬ 
more,  surface  hydroxyl  groups  are  ubiquitous  in  nearly  any 
environment,  as  they  decorate  the  surface  of  metals,  metal 
oxides,  and  organic  materials.  Therefore,  the  primary  focus  of 
the  work  described  below  is  to  provide  fundamental  insight  into 
the  strength  and  structure  of  HD— surface  hydrogen-bonding 
interactions  through  the  systematic  study  of  key  mimics  of  the 
actual  chemical  warfare  agent. 
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Research  on  CWAs  typically  employs  simulants,  which  are 
intended  to  mimic  the  chemical  properties  of  the  live  agent 
without  the  risk  of  exposing  laboratory  personnel  to  extremely 
toxic  compounds.  The  most  common  simulant  for  HD  is  the 
commercially  available  compound  2-CEES.  HD  and  2-CEES 
are  structurally  nearly  identical  and  differ  only  by  an  additional 
terminal  chloro  functional  group  in  HD. 

Hydrogen-bonding  interactions  between  2-CEES  and  the 
free  hydroxyl  groups  on  the  surface  of  a  mixed  Ti02— Si02 
composite  material  have  been  previously  explored  through 
infrared  spectroscopy.4  7  In  those  studies,  the  IR  absorptivity 
of  isolated  SiO— H  stretches  on  the  surface  was  monitored 
during  2-CEES  adsorption.  The  data  revealed  that  the  SiO— H 
band  underwent  extensive  changes,  indicative  of  hydrogen- 
bond  formation  with  both  the  Cl  and  S  moieties  of  the 
adsorbate.  The  three  key  changes  in  the  SiO— H  band  included 
a  shift  to  lower  frequencies,  a  broadening  of  the  band,  and  an 
increase  in  the  absorptivity  of  the  band.  Interestingly,  the 
previous  work  demonstrated  that  the  magnitude  of  each  effect 
was  significantly  different  for  hydrogen  bonding  through  the  Cl 
or  the  S  groups.  This  effect  was  ascribed  to  differences  in  the 
electronic  structure  of  the  moieties,  the  result  of  differing 
degrees  of  hybridization  in  the  lone-pair  electrons  of  the 
hydrogen-bond  acceptor.  One  might  then  speculate  that  these 
differences  may  also  affect  the  strength  of  the  hydrogen  bonds. 

Although  the  hydrogen-bonding  energy  between  a  2-CEES 
or  HD  molecule  and  hydroxyl-containing  compounds  or 
surfaces  has  not  yet  been  (to  the  best  of  our  knowledge) 
reported  in  the  literature,  several  recent  studies  have  focused  on 
the  energetic  and  mechanistic  details  of  hydrogen  bonding  for 
similar  compounds  in  the  gas  phase.8,9  The  understanding  that 
has  emerged  from  those  studies  is  that  the  character  of  a 
hydrogen  bond,  including  its  strength,  bond  length,  and 
spectroscopic  signatures,  is  highly  dependent  on  both  electro¬ 
static  interactions  and  the  extent  of  charge  transfer  between  a 
nonbonding  orbital  of  the  acceptor  and  an  antibonding  orbital 
of  the  hydride  donor.9  For  example,  a  comprehensive 
computational  study8  of  hydrogen-bonding  interactions  at  the 
coB97X-D/aug-cc-pVTZ  level  of  theory  showed  that  the  OH 
stretching  frequency  for  the  hydrogen-bond  donor  in  H20 
shifts  by  88  cm-1  when  it  forms  a  bond  with  the  chloro  group 
of  methyl  chloride.  However,  the  OH  stretching  frequency 
shifts  by  more  than  200  cm-1  when  water  forms  a  hydrogen 
bond  with  the  S  moiety  of  dimethyl  sulfide.  Despite  the  large 
differences  in  the  OH  stretching  frequency  for  these  two 
hydrogen  bonds,  the  calculated  energies  of  the  bonds  differ  by 
only  2  kcal/mol.8  These  computational  investigations  suggest 
that  hydrogen-bond  formation  through  both  the  Cl  and  S 
groups  of  2-CEES  (and  HD)  may  contribute  significantly  to  the 
energetics  of  adsorption  on  hydroxyl-containing  surfaces. 

The  work  presented  here  investigates  the  energy  and 
mechanism  of  hydrogen-bond  formation  between  2-CEES,  as 
well  as  a  series  of  key  test  molecules,  and  the  surface  of  silica. 
Silica  was  chosen  for  this  work  not  only  because  it  is  one  of  the 
most  abundant  materials  found  in  environmental  and  industrial 
settings  but  also  because  it  can  serve  as  an  effective  model 
substrate  with  well-characterized  interfacial  hydroxyl  groups 
that  can  be  prepared  in  a  highly  reproducible  manner.  In  these 
studies,  transmission  infrared  spectroscopic  measurements  are 
used  to  probe  changes  to  the  vibrational  frequency  of  SiO— H 
stretches  during  gas  uptake  and  hydrogen-bond  formation. 
Shifts  in  this  band  provide  insight  into  the  extent  of  charge 
transfer  into  the  antibonding  orbital  of  the  hydroxyl  group 


upon  the  formation  of  the  hydrogen  bond.  However,  the 
binding  energy  depends  on  a  complex  interplay  between  charge 
transfer,  electrostatic  forces,  and  other  intermolecular  forces 
that  cannot  be  predicted  from  spectroscopic  measurements 
alone.  Therefore,  we  employ  temperature-programmed  desorp¬ 
tion  (TPD)  methods  to  interrogate  the  activation  energy  for 
bond  breaking  and  desorption,  which  provides  insight  into  the 
strength  of  the  hydrogen  bond.  Performing  these  studies  under 
ultrahigh  vacuum  (UHV)  conditions  ensures  that  the  silica 
surface  remains  free  from  contamination  throughout  the 
experiment.  Overall,  these  studies  provide  fundamental  insight 
into  hydrogen  bonding  and  the  interactions  at  the  gas— surface 
interface  that  are  responsible  for  HD  uptake,  residence  time, 
and,  ultimately,  surface  reactivity. 

2.  EXPERIMENTAL  SECTION 

Instrumentation.  The  instrumentation  and  general  proce¬ 
dures  used  for  this  project  are  based  on  work  previously 
reported  by  Wilmsmeyer  et  al. 10,11  Experiments  were 
performed  in  an  UHV  chamber  with  a  base  pressure  of  10-9 
Torr.  These  low  pressures  are  necessary  to  minimize  the 
possibility  of  surface  contamination  on  the  silica  sample,  which 
is  a  highly  effective  sorbent  especially  at  cryogenic  temper¬ 
atures.  Control  experiments  were  performed  to  ensure  that 
background  water  and  other  common  UHV  contaminants  did 
not  adsorb  to  the  surface  over  the  duration  of  a  typical 
experiment.  Specifically,  infrared  spectroscopic  measurements 
showed  no  sign  of  changes  to  the  silica  sample  (either  in 
spectral  regions  assigned  to  OH  or  CH  stretches)  over  several 
hours  in  a  vacuum  with  the  surface  at  a  temperature  of  130  K. 
The  total  coverage  of  contaminants  in  this  work  remains  low, 
largely  due  to  the  very  high  surface  area  of  our  particles,  which 
necessitates  very  long  exposure  times  to  accumulate  even  trace 
levels  of  contamination  in  an  UHV  environment. 

Silica  samples  were  prepared  in  a  hydraulic  press  by  packing 
approximately  5  mg  of  Aerosil  (200  m2/g  surface  area,  12  nm 
mean  particle  diameter)  into  a  50  /(m  thick  tungsten  mesh  grid 
(Tech  Etch).  The  mesh  was  mounted  to  a  McAllister  X-Y-Z 
manipulator  via  nickel  clamps  and  5  mm  thick  copper  leads 
connected  to  an  external  power  supply  to  achieve  precise 
control  of  sample  position  and  temperature.  Sample  temper¬ 
ature  was  monitored  with  a  type-K  thermocouple  spot-welded 
to  the  mesh  immediately  adjacent  to  the  silica  sample.  The 
mesh  was  heated  resistively  using  a  custom  external  power 
supply.  The  power  leads  and  thermocouple  wire  pass  through  a 
tube  in  the  manipulator,  which  could  be  filled  with  liquid 
nitrogen  to  cool  the  silica  sample  to  cryogenic  temperatures. 
Prior  to  each  experiment,  the  silica  sample  was  pretreated  by 
heating  to  700  K  for  5  min,  which  removes  possible  trace 
surface  contamination.  Pretreatment  also  significantly  dehy- 
droxylates  the  surface  such  that  it  is  dominated  by  isolated  OH 
groups  at  an  approximate  density  of  ~2  OH/nm2.1"  16  After 
heating,  the  sample  was  cooled  to  cryogenic  temperatures  to 
maximize  uptake  during  dosing. 

Dosing.  2-CEES,  diethyl  sulfide,  chloroethane,  1-chloropro- 
pane,  1-chlorobutane,  and  1-chloropentane  (molecular  struc¬ 
tures  are  provided  in  Figure  l)  were  purchased  from  Sigma- 
Aldrich  and  transferred  to  stainless  steel  cylinders  mounted  on 
a  gas-handling  manifold  configured  with  VCR  seals  and  fittings. 
Each  chemical  was  purified  with  two  or  more  freeze— pump- 
thaw  cycles  immediately  prior  to  use.  All  compounds  employed 
for  this  work  are  liquids  at  room  temperature  but  have 
sufficiently  high  vapor  pressures  to  deliver  gas-phase  molecules 


366 


D0!:10.1021/jp509516x 
J.  Phys.  Chem.  C2015,  119,  365-372 


Article 


The  Journal  of  Physical  Chemistry  C 


Cl 


HD 


Cl 


2-CEES 


cr  v  v 

Chloropentane 


cr  v  ^ 

Chlorobutane 


Cl 


form  surface  coverage  is  evidenced  by  poor  alignment  of  the 
trailing,  high-temperature  edges  of  the  TPD  distributions.17 

Temperature-Programmed  Desorption.  TPD  experi¬ 
ments  were  performed  by  heating  the  sample  and  detecting 
desorbed  species  with  a  doubly  differentially  pumped,  axially 
mounted  quadrupole  mass  spectrometer  (Extrel)  tuned  to  the 
most  abundant  mass  fragments  for  each  molecule.  The  mass 
spectrometer  views  a  0.35  cm2  spot  on  the  sample  with  an 
acceptance  angle  of  ±0.4°  such  that  the  signal  from  molecules 
that  may  desorb  from  the  sample  holder  is  minimized.  The 
heating  rate  for  the  TPD  measurements  was  maintained  at  0.2 
K/s  with  a  proportional-integral-derivative  controller  (Honey¬ 
well)  and  custom  power  supply. 


Diethyl  Sulfide 
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3.  RESULTS 

Infrared  Studies  of  Adsorption.  Previous  studies  of  2- 
CEES  and  other  polar  compounds  demonstrated  the  use  of 
infrared  spectroscopy  in  the  characterization  of  interactions 
between  silica  and  adsorbed  molecules.4,8  Figure  2  shows  a 


Chloroethane 

Figure  1.  Schematics  for  HD  and  molecules  explored  in  this  work. 


to  the  UHV  chamber  through  a  directional  capillary  array  doser 
positioned  approximately  5  mm  from  the  silica  surface.  The 
doser  was  mounted  on  a  translation  stage  and  retracted  when 
not  in  use.  The  manifold  was  kept  under  vacuum  at  all  times 
except  during  dosing  to  minimize  contamination  and  was 
vacuum-heated  for  12  h  at  over  500  K  to  clean  the  system 
between  experiments  with  each  molecule.  The  temperature  of 
the  silica  sample  during  dosing  was  adjusted  for  each  molecule 
to  attain  high  (at  least  75%)  surface  coverage  within  3—7  min. 

Infrared  Spectroscopy.  The  silica  sample  was  monitored 
with  transmission  Fourier- transform  infrared  spectroscopy 
(FTIR)  before,  during,  and  after  dosing  to  identify  adsorbed 
molecules  and  changes  in  surface  vibrational  modes.  The  IR 
spectra  were  recorded  using  a  Nicolet  Nexus  670  FTIR 
spectrometer  integrated  with  the  UHV  chamber.  All  IR  spectra 
consist  of  128  scans  recorded  at  4  cm-1  resolution.  IR  spectra 
of  clean  silica  employed  an  empty  spot  on  the  mesh  as  the 
reference  background.  Spectra  of  species  adsorbed  on  the  silica 
sample  employed  a  spectrum  of  the  clean  silica  as  a  reference. 

Annealing.  Previous  work  indicates  that  effective  TPD 
measurements  require  a  homogeneous  distribution  of  adsorbed 
molecules  throughout  the  particulate  sample.17  We  achieved 
this  distribution  by  annealing  the  silica  sample  immediately 
after  dosing  to  a  sufficiently  high  temperature  for  the  adsorbed 
molecules  to  diffuse  evenly  through  the  packed  bed  of  particles. 
A  different  anneal  temperature  was  chosen  for  each  molecule  to 
reduce  the  surface  concentration  to  submonolayer  coverage  and 
achieve  uniform  distributions.  Anneal  temperatures  were  219  K 
for  2-CEES,  198  K  for  diethyl  sulfide,  163  K  for  chloropropane, 
181  K  for  chlorobutane,  and  193  K  for  chloropentane. 
Chloroethane  did  not  require  annealing  to  remove  multilayers 
and  achieve  a  uniform  distribution  on  the  surface.  Surface 
coverage  for  all  molecules  was  controlled  by  varying  the 
annealing  duration.  After  annealing,  the  samples  were  rapidly 
cooled  to  well  below  their  desorption  temperatures,  where  the 
TPD  measurements  were  initiated.  Sufficient  annealing  can  be 
confirmed  experimentally  by  comparing  plots  of  desorption  rate 
vs  temperature  at  different  experimental  coverages.  Nonuni- 


Figure  2.  Infrared  spectrum  of  a  clean  silica  sample  following 
pretreatment  at  700  K  for  5  min  (see  the  text). 


spectrum  of  clean  silica  immediately  before  exposure  to  the 
adsorbate.  The  sharp  absorption  feature  at  3750  cm-1  is 
characteristic  of  free  surface  silanol  groups.18  The  shoulder  at 
3700  cm-1  is  indicative  of  geminal  silanols,  but  the  small  size  of 
this  band  relative  to  the  main  jaeak  indicates  that  isolated 
silanols  dominate  the  surface.12-1' ,19  This  is  consistent  with  a 
study  that  measured  the  distributions  of  free,  geminal,  and 
hydrogen-bonded  OH  surface  groups  per  unit  area  as  a 
function  of  temperature  from  473  to  1473  K.16  The 
wavenumber  region  from  1400  to  650  cm-1  is  dominated  by 
bulk  infrared  modes  of  silica  that  obscure  several  adsorbate 
bands  (see  below),  including  those  associated  with  the  scissor 
mode  for  CH2  adjacent  to  Cl  (1215  cm-1)  and  the  C— Cl 
stretch  (600—800  cm  :)  for  some  of  the  molecules  used  in  this 
work. 

IR  difference  spectra  following  uptake  of  chlorobutane, 
diethyl  sulfide,  and  2-CEES  on  silica  are  presented  in  Figure  3A. 
Figure  3B  shows  IR  difference  spectra  for  the  silica  particles 
following  uptake  of  the  aliphatic  chloroalkanes.  The  inverted 
feature  at  3750  cm-1  is  the  result  of  a  shift  in  that  band  as  the 
surface  silanol  groups  transform  from  free  SiOH  to  hydrogen- 
bonded  groups.  With  the  assumption  that  the  absorbance 
intensity  of  the  free  OH  band  is  directly  proportional  to  the 
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Figure  3.  (A)  Infrared  spectra  of  adsorbed  diethyl  sulfide  (top,  red),  chlorobutane  (middle,  green),  and  2-CEES  (bottom,  blue)  on  silica.  The  spectra 
are  normalized  to  the  integrated  area  of  the  free  OH  band  such  that  the  spectra  represent  equivalent  effective  SiOH  occupation  fractions,  9  (see  the 
text).  (B)  Infrared  spectra  of  adsorbed  chloroalkanes. 


surface  concentration  of  free  OH,  the  decrease  in  intensity  of 
the  feature  at  3750  cm-1  following  exposure  to  the  gas  of 
interest  is  used  to  determine  the  fraction  of  silanol  groups 
involved  in  hydrogen  bonding.  In  this  study,  surface  coverage, 
6,  is  reported  in  terms  of  the  fraction  of  silanol  groups  bound  to 
an  adsorbate;  for  example,  6=1  represents  occupation  of  all 
hydrogen-bonding  sites.  The  monotonic  profile  and  excellent 
alignment  of  the  TPD  curves  shown  in  Figure  4  are  evidence 


Temperature  (K) 


Figure  4.  Experimental  (blue  circles)  and  simulated  (black  lines)  TPD 
curves  for  2-CEES  at  different  coverages.  From  top  to  bottom:  9  = 
0.79,  0.46,  0.32,  0.21,  and  0.17.  The  inset  shows  a  6th-order 
polynomial  fit  to  ^  error  values  between  simulated  and  experimental 
data  used  to  determine  the  best  prefactor  value. 


for  the  absence  of  significant  multilayers.  Furthermore,  the 
decrease  in  IR  signal  from  hydrogen-bonded  silanol  groups 
during  desorption  directly  correlates  with  a  decrease  in 
adsorbate  IR  bands  and  a  rise  in  mass  spectrometer  (MS) 
signal  at  the  parent  molecular  mass.  This  correlation  is  evidence 
that  the  adsorbed  molecules  are  bound  through  the  SiOH 
groups. 

The  large  shift  in  the  silanol  absorption  band,  clearly  evident 
in  Figures  3  and  5,  indicates  that  hydrogen  bonds  are  an 
important  driving  force  in  the  uptake  of  these  molecules.  Three 
key  characteristics  of  the  IR  band  for  the  SiO— H  stretch  stand 
out  in  the  spectra  for  adsorbed  molecules:  (l)  the  SiO— H 
stretch  band  shifts  to  lower  wavenumber  (red-shift),  (2)  the 
hydrogen-bonded  OH  band  has  a  much  larger  intensity  than 
the  free  OH  band,  and  (3)  this  band  broadens  significantly 


4000  3000  2000  1000 

Wavenumber  (cm1) 

Figure  S.  Infrared  spectra  of  gas-phase  2-CEES  (black)  and  2-CEES 
adsorbed  on  silica  (blue).  The  inset  displays  IR  bands  corresponding 
to  methyl  and  methylene  stretching  modes. 


upon  adsorption.  The  red-shift  is  due  to  charge  transfer  as  the 
hydrogen-bond  acceptor  donates  electron  density  to  SiOH. 
Upon  hydrogen-bond  formation,  n  electrons  from  the  hydro- 
gen-bond  acceptor  interact  with  the  a*  antibonding  orbital  on 
the  silanol,  thereby  weakening  the  O— H  bond  and  adding 
significant  anharmonicity  to  the  stretching  mode.20,21  The 
increase  in  peak  intensity  is  due  to  the  greater  oscillator 
strength  of  hydrogen-bonded  donor  stretches.22  Finally,  the 
band  broadening  occurs  because  of  the  inhomogeneity  in  the 
configuration  of  the  hydrogen-bonded  molecules,  which  affects 
the  strength  and  anharmonicity  of  each  bond  in  different  ways. 
Furthermore,  coupling  to  other  modes  may  reduce  the  lifetime 
of  the  O— H  vibration,  which  increases  the  homogeneous  band 
broadening.23-26 

The  bands  near  3000  cm-1  are  assigned  to  symmetric  and 
asymmetric  alkane  stretching  modes  corresponding  to  mole¬ 
cules  on  the  surface.  A  comparison  of  methylene  stretching 
modes  for  each  adsorbed  molecule  to  corresponding  gas-phase 
modes  for  the  same  molecule  provides  insight  into  the 
adsorption  mechanism.  Figure  5  shows  spectra  of  2-CEES 
adsorbed  on  silica  and  in  the  gas  phase.  The  relative  integrated 
areas  of  each  band  are  similar,  which  indicates  molecular 
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Figure  6.  (A)  Ed(0)  curves  for  the  desorption  of  2-CEES,  diethyl  sulfide,  and  chlorohutane  from  partially  hydroxylated  silica.  (B)  Ed(0)  curves  for 
the  desorption  of  four  chloroalkanes  from  partially  hydroxylated  silica. 


adsorption  without  subsequent  decomposition  or  reaction. 
However,  the  bands  assigned  to  t/a(CH2)  and  trs(CH3)  shift 
upon  adsorption  by  5  and  10  cm-1,  respectively,  compared  to 
the  same  modes  in  the  gas  phase.  Similar  shifts  are  observed  in 
all  molecules  investigated  here.  They  are  likely  due  to 
dispersion  interactions  with  the  surface,  restricted  molecular 
degrees  of  freedom  for  the  adsorbates,  and  changes  in  the  CH2 
bending  modes  (the  first  overtone  is  in  Fermi  resonance  with 
the  CH2  stretches).  Unfortunately,  the  CH2  bending  modes 
could  not  be  detected  in  the  present  work  due  to  obstruction 
by  the  high  absorbance  of  silica. 

Temperature-Programmed  Desorption.  The  desorption 
energies  were  determined  by  using  TPD.  After  dosing  and 
annealing,  the  silica  sample  was  heated  at  a  constant  rate  of  0.2 
K/s.  This  rate  was  chosen  to  allow  adsorbed  molecules 
sufficient  time  to  diffuse  throughout  the  pressed  sample  before 
desorbing.  The  mass  spectrometer,  tuned  to  mass  fragments 
unique  to  the  desorbing  molecule,  tracked  changes  in  the  gas- 
phase  number  density  of  the  analytes  as  a  function  of  surface 
temperature.  The  raw  signal  was  then  multiplied  by  T1/2  to 
correct  for  the  velocity-dependent  detection  probability  that 
affects  line-of-sight,  differentially  pumped,  mass  spectrometric 
TPD  measurements.27  The  resulting  distribution  was  directly 
proportional  to  desorption  rate.  Molecular  desorption,  as 
opposed  to  decomposition  on  the  surface  and  subsequent 
desorption  of  products,  was  confirmed  by  comparing  the 
fragmentation  pattern  for  the  pure  gases  to  that  recorded 
during  desorption.  The  desorption  rates  for  2-CEES  at  different 
initial  coverages  are  shown  in  Figure  4. 

The  variation  of  the  desorption  rate  with  the  initial  coverage 
provides  a  qualitative  understanding  of  bonding.  The  alignment 
of  the  trailing  edges  indicates  that  adsorbed  molecules  have 
sufficient  mobility  to  diffuse  to  the  strongest  available 
adsorption  sites  during  the  temperature  ramp  prior  to 
desorption.17  The  maximum  in  the  desorption  rate  shifts  to 
higher  temperatures  with  lower  coverage.  This  is  evidence  that 
the  residence  time  and  binding  is  governed  by  molecule- 
surface,  rather  than  molecule— molecule,  forces. 

Quantitative  interpretation  of  TPD  spectra  is  based  on  the 
Polanyi— Wigner  relationship  (eq  l),  which  can  be  rewritten  as 
eq  2  to  express  desorption  energy  as  a  function  of  coverage. 
The  experimental  TPD  data  were  analyzed  by  using  the 
inverted  Polanyi— Wigner  equation  to  determine  the  desorption 


energy  at  each  temperature  and  corresponding  coverage 
according  to  the  methods  outlined  by  others17  and  recently 
used  to  study  the  desorption  energetics  of  molecules  from 
particulate  silica.10'11  The  desorption  energies  are  then  plotted 
as  a  function  of  coverage  and  modeled  with  first-order  kinetics. 
Numerical  integration  of  the  desorption  energy  ( Ed )  vs 
coverage  (#)  curve  obtained  from  the  highest  experimental 
coverage  gives  simulated  TPD  curves  for  any  desired  initial 
coverage  lower  than  that  highest  experimental  coverage.  The 
prefactor  ( v  in  eq  l)  is  assumed  to  be  independent  of  coverage 
and  temperature  for  this  type  of  analysis  and  is  used  as  a  fitting 
parameter  to  maximize  agreement  between  simulated  and 
experimental  desorption  rate  curves.  Specifically,  the  prefactor 
is  adjusted  to  minimize  the  sum  of  the  squared  residuals 
between  experimental  data  and  simulated  data  at  each  coverage. 
Figure  4  shows  simulated  data  (black  lines)  obtained  with  an 
optimized  prefactor  plotted  over  the  experimental  Ed  vs 
temperature  data  for  five  initial  coverages.  The  inset  shows 
the  SSR  plot  used  to  optimize  the  prefactor. 

Ts)  =  i>(0,  Ts)  exp [— Ed(0)/fcBTs]0" 
at  w 

Ea(S)  =  (2) 

After  optimizing  the  Polanyi— Wigner  prefactor,  the  reported 
desorption  energy  is  determined  from  the  plot  of  desorption 
energy  vs  coverage  (Figure  6A).  The  curvature  in  this  plot  is 
likely  due  to  a  nonuniformity  of  the  silica  surface.  Some  silanol 
sites,  defect  sites  in  particular,  likely  contribute  higher  binding 
energy  sites  to  the  surface.  Furthermore,  hydrogen  bonding, 
though  the  most  significant  force  in  the  physisorption  of  these 
molecules,  is  not  the  only  type  of  interaction  that  determines 
the  overall  binding  energy.  Other  surface— adsorbate  inter¬ 
actions,  including  dispersion  forces,  can  significantly  affect  the 
activation  energy  for  desorption.  Thus,  the  local  variations  in 
silanol  density  likely  affect  the  desorption  energy  for  individual 
molecules.  In  this  study,  our  focus  is  on  the  desorption  energy 
of  a  single  adsorbate  on  a  defect-free  region  of  the  silica  surface; 
therefore,  the  linear  portion  of  the  Ed  curve  is  extrapolated  to 
zero  coverage.  Desorption  energy  distributions  for  2-CEES, 
diethyl  sulfide,  and  chlorobutane  are  presented  in  Figure  6A. 
Figure  6B  contains  desorption  energy  distributions  for  the 
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chloroalkanes.  The  extrapolated,  zero-coverage  limit  desorption 
energies  for  each  of  the  molecules  studied  are  presented  in 
Table  1. 


Table  1.  Desorption  Energies  and  Pre-Exponential  Factors 
for  Molecules  Tested  in  This  Study” 

molecule 

Ed  (±2.3  kj/mol) 

logM  (±0.6) 

2-CEES 

42.9 

7.7 

diethyl  sulfide 

37.5 

6.5 

chloropentane 

43.7 

9.1 

chlorobutane 

35.4 

7.3 

chloropropane 

29.2 

6.3 

chloroethane 

24.4 

5.9 

aError  values  are  ±  one  standard  deviation 
independent  TPD  experiments. 

calculated  from  three 

4.  DISCUSSION 

We  have  explored  the  fundamental  nature  of  gas— surface 
hydrogen-bond  formation  between  simulant  molecules  (2- 
CEES,  diethyl  sulfide,  and  chlorobutane)  and  partially 
hydroxylated  silica  to  help  construct  an  understanding  of  how 
sulfur  mustard  gas  adheres  to  hydroxyl-containing  surfaces.  In 
this  study,  infrared  spectroscopy  was  employed  to  probe  the 
extent  of  charge  transfer,  a  key  mechanistic  step  in  hydrogen- 
bond  formation,  from  the  adsorbate  to  the  hydroxyl  group  of  a 
well-characterized  silica  surface.  The  overall  strengths  of  the 
hydrogen  bonds  were  gauged  by  performing  TPD  measure¬ 
ments  to  determine  the  activation  energy  for  desorption. 
Together,  these  techniques  provide  insight  into  the  roles  that 
dipole— dipole  interactions,  charge  transfer,  and  dispersion 
forces  play  in  the  uptake  of  mustard  on  polar  surfaces. 

2-CEES-Silica  Hydrogen-Bond  Formation.  2-CEES 
uptake  on  silica  is  primarily  driven  by  the  formation  of 
hydrogen  bonds  with  surface  hydroxyl  groups.  These  hydrogen¬ 
bonding  interactions  are  evidenced  in  the  IR  spectra  of 
adsorbed  2-CEES  (Figure  5),  which  clearly  show  the  reduction 
of  the  band  associated  with  free  surface  hydroxyl  groups 
(reduced  absorbance  at  3750  cm-1)  as  two  new  broad  IR 
features,  signatures  of  hydrogen  bonding,  emerge.  The  two 
broad  but  distinct  IR  absorbance  features  at  3598  and  3359 
cm-1  are  in  the  infrared  spectral  range  associated  with  surface 
SiO— H  stretches  but  are  red-shifted  significantly  from  the  free 
silanol  band.  This  red-shift  reflects  a  decrease  in  the  SiO— H 
vibrational  energy  caused  by  charge  transfer  from  the  hydrogen- 
bond  acceptor  to  the  antibonding  orbital  of  the  SiO— H 
bond.8,28  In  this  way,  infrared  spectroscopic  measurement  of 
the  SiO— H  vibrational  frequency  is  a  probe  of  the  extent  of 
charge  transfer  in  gas— surface  hydrogen-bond  formation  and 
reveals  two  distinct  types  of  charge  donors  (lone-pair 
hydrogen-bond  acceptors)  within  2-CEES. 

The  chlorine  and  sulfur  moieties  within  the  2-CEES 
molecule  are  both  potential  hydrogen-bond  acceptor  sites. 
Therefore,  we  have  studied  the  uptake  of  chlorobutane  and 
diethyl  sulfide  to  help  identify  the  contributions  from  each  site. 
Spectra  for  adsorbed  chlorobutane  show  a  strong  absorbance 
feature  at  3571  cm-1,  which  is  due  to  excitation  of  the  SiO— H 
stretch  associated  with  the  C1---H— OSi  hydrogen  bond.  Diethyl 
sulfide  also  forms  hydrogen  bonds  with  the  surface,  but  the 
band  associated  with  the  SiO— H  stretch  is  much  further  red- 
shifted  in  the  IR  spectra.  This  band  appears  at  3280  cm-1. 
Thus,  we  assign  the  3598  and  3359  cm-1  absorption  features  in 


the  IR  spectra  of  adsorbed  2-CEES  to  Cl-H— OSi  and  S---H— 
OSi  hydrogen  bonds,  respectively,  in  agreement  with  a  previous 
study.4 

The  infrared  spectral  assignments  reveal  that  the  extent  of 
charge  transfer  during  uptake  and  hydrogen-bond  formation  is 
greater  for  molecules  that  bind  through  the  sulfur  than  through 
the  chlorine  constituents  of  the  molecule  and  that  both  types  of 
interactions  occur  for  2-CEES.  However,  the  extent  of  charge 
transfer  does  not  always  correlate  with  the  strength  of  the 
hydrogen  bond.  Electrostatic  interactions  also  play  a  major  role 
in  hydrogen-bond  formation,  and  further  measurements  are 
necessary  to  assess  the  relative  contribution  of  each  type  of 
hydrogen  bond  to  the  overall  adsorption  process.  Therefore,  we 
have  performed  a  series  of  TPD  studies  to  help  determine  the 
overall  energy  required  to  rupture  the  hydrogen  bonds  and 
drive  the  molecules  from  the  surface. 

Charge  Transfer  and  Electrostatic  Interactions.  Pre¬ 
vious  studies  in  solution,  in  the  gas  phase,  and  at  the  gas- 
surface  interface  reveal  a  positive,  linear  correlation  between  At/ 
(the  shift  in  the  frequency  for  the  O— H  stretch  of  a  hydrogen 
bond)  and  intermolecular  hydrogen-bond  strength  for 
electronically  similar  moieties. 10,11,29-33  On  the  basis  of  this 
previous  work  alone,  one  would  infer  that  the  larger  Au  in 
spectra  of  adsorbed  diethyl  sulfide  (At/  =  485  cm-1)  compared 
to  chlorobutane  ( Au  =  190  cm-1)  is  indicative  of  a  stronger 
hydrogen  bond.  However,  as  described  above,  electrostatic  (in 
this  case,  dipole— dipole)  interactions  also  play  a  role  in 
hydrogen-bond  forces.  Calculations  indicate  that  the  molecular 
dipole  moments  for  chlorobutane  and  diethyl  sulfide  are  2.47 
and  1.68  D,  respectively.34  Thus,  on  the  basis  of  the  magnitude 
of  the  dipoles,  one  might  predict  that  chlorobutane  would  form 
a  stronger  hydrogen  bond  with  silanol  than  diethyl  sulfide.  This 
is  contrary  to  the  conclusion  based  solely  on  the  extent  of  the 
red-shift  in  the  SiO— H  vibrational  frequency  in  the  hydrogen- 
bonded  complexes. 

The  interplay  between  dipole— dipole  interactions  and  charge 
transfer  has  been  explored  in  a  recent  computational  study  on 
hydrogen  bonding  of  dimethyl  sulfide  and  chloromethane  to 
gas-phase  water  (among  a  variety  of  other  species).8  In  that 
study,  Au  for  the  O— H  stretch  was  125  cm-1  larger  for  water 
bound  to  dimethyl  sulfide  than  for  water  bound  to  chloro¬ 
methane.  However,  the  calculated  hydrogen-bond  strengths  for 
both  molecules  were  similar  relative  to  other  molecules  studied 
in  that  work.  Analysis  of  our  TPD  measurements  for  diethyl 
sulfide  and  chlorobutane  (See  Figure  6A  and  Table  l)  is 
consistent  with  these  calculations.  That  is,  the  activation 
energies  for  desorption  for  diethyl  sulfide  and  chlorobutane 
from  silica  were  found  to  be  nearly  the  same  within  the  error 
bars  for  these  measurements  (37.5  ±  2.3  and  35.4  ±  2.3  kj/ 
mol).  Thus,  we  conclude  that  differences  in  the  charge  transfer 
contribution  are  balanced  by  other  intermolecular  forces,  such 
as  dipole— dipole  interactions,  to  yield  similar  activation 
energies  for  desorption.  Therefore,  in  the  case  of  2-CEES 
adsorption,  hydrogen-bond  formation  through  the  sulfur  or  the 
chlorine  component  likely  leads  to  very  similar  binding 
strengths  (similar  to  those  of  diethyl  sulfide  and  chlorobutane) 
for  isolated  1:1  type  interactions.  However,  if  2-CEES  were  to 
interact  simultaneously  with  an  SiOH  group  through  the 
chlorine  and  another  SiOH  group  through  the  sulfur,  then  the 
activation  energy  for  desorption  may  be  much  higher. 

2-CEES  Adsorption  Mechanism.  We  have  investigated 
whether  each  adsorbed  2-CEES  molecule  interacts  with  a  single 
silanol  group  (1:1  adsorbate:silanol  ratio)  or  two  adjacent 
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silanol  groups  (1:2  ratio).  As  described  above,  the  surface 
hydroxyl  density  on  silica  pretreated  to  700  K  is  estimated  to  be 
~2  OH/nm2,  which  corresponds  to  an  approximate  spacing  of 
0.7  nm  between  silanols.3,16'35  This  spacing  is  likely  sufficient  to 
allow  the  two  polar  functional  groups  in  2-CEES,  separated  by 
an  ethyl  chain,  to  simultaneously  bind  to  adjacent  hydroxyls. 
However,  analysis  of  the  relative  surface  coverages  for  2-CEES 
and  diethyl  sulfide  at  identical  SiOH  occupation  densities 
suggests  that  2-CEES  adsorbs  similarly  to  diethyl  sulfide;  that  is, 
there  is  only  one  molecule— surface  hydrogen  bond  per 
adsorbed  molecule.  More  specifically,  the  spectra  presented  in 
Figure  3A  are  normalized  such  that  they  represent  the  same 
number  of  occupied  SiOH  surface  sites.  At  this  coverage,  we 
find  that  the  integrated  areas  under  the  bands  for  the  CH2 
asymmetric  stretching  mode  of  each  molecule  are  nearly 
identical,  strongly  suggesting  that  the  adsorbate-to-SiOH  ratios 
(necessarily  1:1  for  diethyl  sulfide)  are  the  same.  This 
conclusion  is  supported  by  the  TPD  data,  which  indicates 
that  the  desorption  energy  for  2-CEES  is  only  15%  greater  than 
that  of  diethyl  sulfide,  as  opposed  to  the  much  larger  energy 
increase  one  might  expect  if  two  hydrogen  bonds  were  formed 
per  molecule. 

Stabilization  from  Secondary  Interactions.  The  slightly 
higher  desorption  energy  for  2-CEES  compared  to  those  of 
chlorobutane  or  diethyl  sulfide  is  likely  due  to  surface- 
adsorbate  interactions  other  than  the  hydrogen  bonds.  Weaker 
secondary  interactions  likely  play  a  role  in  controlling  the 
activation  energy  for  desorption.  Previous  work  reported  a 
binding  energy  increase  of  7  kj/mol  per  additional  methylene 
for  a  series  of  linear  hydrocarbons  on  MgO(lOO).36  We  have 
observed  a  similar  trend  on  amorphous  silica.  Figure  6A  shows 
the  coverage-dependent  desorption  energies  for  a  series  of 
linear  chloroalkanes  that  differ  only  by  alkane  chain  length.  As 
with  2-CEES,  diethyl  sulfide,  and  chlorobutane,  the  desorption 
energy  is  largely  invariant  with  coverage,  except  for  the  lowest 
coverages  due  to  a  small  concentration  of  sites  that  lead  to 
unusually  high  binding  energies.  The  linear  portions  of  the 
Ed(0)  distributions  are  extrapolated  to  the  zero  coverage  limit 
and  the  values  are  provided  in  Table  1.  The  chloroalkane 
desorption  energies  range  from  24  kj/mol  for  the  shortest 
chain  molecule  studied  to  44  kj/mol  for  the  longest  chain 
molecule.  Each  methylene  group  contributes  5—8  kj/mol  to  the 
overall  desorption  energy.  It  is  unlikely  that  this  energy 
difference  is  caused  by  changes  to  the  SiO— H---C1  hydrogen- 
bond  strength,  because  the  fundamental  nature  of  the  hydrogen 
bonds,  as  gauged  by  the  SiO— H  stretching  frequency  (A v  = 
180  cm-1  for  all  four  molecules),  appears  to  be  identical  (see 
Figure  3B).  Furthermore,  calculations  indicate  that  the 
magnitude  of  the  dipoles  for  the  chloroalkanes  are  not  strongly 
dependent  on  chain  length.  The  permanent  dipole  moments 
range  from  2.3  D  for  chloroethane  to  2.5  D  for 
chloropentane.34  Thus,  the  measured  increase  in  desorption 
energy  as  a  function  of  chain  length  likely  arises  from  surface- 
adsorbate  dispersion  interactions.  On  the  basis  of  these  results, 
we  reason  that  similar  dispersion  interactions  for  the  larger  2- 
CEES  molecule,  relative  to  diethyl  sulfide  and  chlorobutane, 
account  for  the  differences  in  the  desorption  energies  for  these 
molecules.  Finally,  we  note  that,  in  addition  to  the  desorption 
energy,  the  Polanyi— Wigner  pre-exponential  factors  ( v )  for  this 
series  of  chloroalkanes  depend  on  chain  length.  This 
phenomenon  is  well-known  from  other  work36  and  is  likely 
due  to  the  increased  entropic  penalty  of  constraining  larger 


molecules  to  a  surface  as  well  as  differences  in  the  diffusivity  of 
the  molecules  through  the  silica  particles.36 

5.  SUMMARY 

We  have  found  that  the  uptake  of  2-CEES  on  silica  is  driven 
primarily  by  the  formation  of  hydrogen  bonds  to  the  thioether 
and  chloro  moieties.  Although  both  types  of  hydrogen  bonds 
have  similar  overall  strength,  the  contributions  to  bond  strength 
from  other  forces,  including  charge  transfer,  are  apparently 
different.  The  S---H—  OSi  bond  has  a  larger  charge  transfer 
contribution,  whereas  the  Cl  -H— OSi  bond  has  a  larger 
contribution  from  other  intermolecular  forces.  Dispersion 
forces  between  the  hydrocarbon  segments  of  the  adsorbed 
molecules  and  surface  atoms  provide  additional  stabilization. 
Each  adsorbed  2-CEES  molecule  occupies  a  single  silanol 
group,  even  though  the  sulfur  and  chlorine  atoms  have 
sufficient  spacing  to  span  adjacent  silanol  groups.  The  1:1 
bonding  motif  may  be  a  result  of  a  significant  entropic  penalty 
that  would  be  required  for  dual  hydrogen-bond  formation. 

A  major  motivation  of  this  work  was  to  develop  insight  into 
the  desorption  energetics  and  binding  mechanisms  for  the 
chemical  warfare  agent  HD  on  the  surface  of  an  environ¬ 
mentally  relevant  hydroxylated  material.  HD  has  the  same 
molecular  structure  as  2-CEES  except  for  the  addition  of  a 
terminal  chlorine  atom;  therefore,  it  is  likely  that  HD  and  2- 
CEES  have  similar  binding  mechanisms  and  desorption 
energies.  The  extra  heteroatom  may  increase  desorption  energy 
somewhat  due  to  added  surface— adsorbate  dispersion  inter¬ 
actions,  but  this  will  be  highly  dependent  on  the  geometry  of 
the  adsorbed  molecules  and  degree  of  silica  hydroxylation. 
Notwithstanding,  we  predict  that  the  overall  desorption  energy 
for  HD  from  silica  will  be  within  several  kilojoules  per  mole  of 
the  values  reported  here  for  its  closest  simulant,  2-CEES. 
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